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Abstract. A few seconds after bounce in a core-collapse supernova, the shock wave 
passes the density region corresponding to resonant neutrino oscillations with the 
"atmospheric" neutrino mass difference. The transient violation of the adiabaticity 
condition manifests itself in an observable modulation of the neutrino signal from a 
future galactic supernova. In addition to the shock wave propagation effects that were 
previously studied, a reverse shock forms when the supersonically expanding neutrino- 
driven wind collides with the slower earlier supernova ejecta. This implies that for 
some period the neutrinos pass two subsequent density discontinuities, giving rise to a 
"double dip" feature in the average neutrino energy as a function of time. We study 
this effect both analytically and numerically and find that it allows one to trace the 
positions of the forward and reverse shocks. We show that the energy dependent 
neutrino conversion probabilities allow one to detect oscillations even if the energy 
spectra of different neutrino flavors are the same as long as the fluxes differ. These 
features are observable in the Ve signal for an inverted and in the v,, signal for a normal 
neutrino mass hierarchy, provided the 13-mixing angle is "large" (sin^ i?i3 ^ 10^^). 



PACS numbers: 14.60.Pq, 97.60.Bw 



1. Introduction 

While galactic supernovae are rare, the proliferation of existing or proposed large 
neutrino detectors has considerably increased the confidence that a high-statistics 
supernova (SN) neutrino signal will eventually be observed. The scientific harvest would 
be immense. Most importantly for particle physics, the detailed features of the neutrino 
signal may reveal the nature of the neutrino mass ordering that is extremely difficult to 
determine experimentally. 
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On the other hand, a detailed measurement of the neutrino signal from a galactic 
SN could yield important clues of the SN explosion mechanism. Neutrinos undoubtedly 
play a crucial role for the SN dynamics, and neutrino energy deposition behind the SN 
shock is able to initiate and power the SN explosion 0I2|- It is, however, still unclear 
whether this energy deposition is indeed sufficiently strong, and current state-of-the- 
art models still have problems to produce robust explosions with the observed energies 
(cf., e.g., Ref. 0). Empirical constraints on the physics deep inside the SN core would 
therefore be extremely useful, and neutrinos are the only way for a direct access besides 
gravitational waves j3]. 

The neutrinos emitted by the collapsed SN core will pass through the mantle 
and envelope of the progenitor star and on the way encounter a vast range of 
matter densities p from nearly nuclear at the neutrinosphere to that of interstellar 
space. The Wolfenstein effect |^ causes a resonance in neutrino oscillations \^ when 
Ami cos2'd /2Ey = ±\/2Gf1^p, where the plus and minus sign refers to neutrinos u and 
antineutrinos u, respectively. Therefore, depending on the sign of Am^, the resonance 
occurs in the u or the u channel [7j. For the "solar" neutrino mass-squared difference 
of Am^i ^ 81 meV^ |H] one refers to the "L-resonance" (low density) while for the 
"atmospheric" one of lAm^gl ~ 2300 meV^ [9j to the "H-resonance" (high density). 
The resonance is particularly important for 13-oscillations because the 13-mixing angle 
is known to be small so that the classic MSW enhancement of flavor conversion in an 
adiabatic density gradient is a crucial feature [Tn| . 

The passage of the shock wave through the density of the H-resonance a few seconds 
after core-bounce (Fig. may break adiabaticity, thereby modifying the spectral 
features of the observable neutrino flux. Therefore, it is conceivable that a future large 
neutrino detector can measure a small modulation of the neutrino signal caused by the 
shock-wave propagation, an effect first discussed by Schirato and Fuller in a seminal 
paper [TTj and elaborated by a number of subsequent authors [121 CHI IHI • Since the 
density of the H-resonance depends on energy, the observation of such a modulation 
in different neutrino energies would allow one to trace the shock propagation. On the 
other hand, the occurrence of this effect depends on the sign of Am|]^ and the value of 
■(^is, so that observing it in the Ue spectra, the experimentally most accessible channel, 
would imply that the neutrino mass ordering is inverted and that sin^'i^is ^ 10~^. 

We here explore a new feature of the shock-wave "fingerprint" in the neutrino 
signal. Some time after the onset of the explosion a neutrino-driven baryonic wind 
develops and collides with the earlier, more slowly expanding supernova ejecta. A 
reverse shock is seen to form in all models which were computed with sufficient resolution 
and seems to be a generic feature, although the exact propagation history depends on 
the detailed dynamics during the early stages of the supernova explosion. Moreover, 
violent convective instabilities and large anisotropies are observed in the neutrino-heated 
layer behind the SN shock in multi-dimensional SN simulations and imply a significant 
angular variation of the instantaneous density profiles even in a single star. Still, the 
simultaneous propagation of a direct and an inverse shock wave imply that for some 
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Figure 1. Shock and reverse-shock propagation in our numerical calculation. The 
density profile is shown at the indicated instances after core bounce. The density 
region pn corresponds to resonant neutrino oscillations with the atmospheric mass 
difference, pl to the solar one. The width of the bands reflects the expected energy 
range of SN neutrinos. 

time the neutrinos may encounter two subsequent density discontinuities, leading to 
significantly different spectral features than are caused by a single crossing. 

Numerical SN simulations do not yet lead to robust explosions so that it remains 
unclear if an important physical ingredient is missing in our understanding of the SN 
phenomenon. Even assuming that only minor aspects need to be tweaked, this situation 
implies that hydrodynamic state-of-the-art models of successful explosions with self- 
consistently determined neutrino fluxes and spectra do not exist at present. Therefore, 
we limit our discussion to a simple analytic study with a schematic density profile 
and an example of a detailed numerical model, using for both cases schematic primary 
neutrino fiuxes and spectra that are representative for neutrino spectra discussed in the 
literature. In this way we identify a generic "double dip" feature caused by the presence 
of two density discontinuities. We find that the position of the two dips in time can 
be connected to the positions of the forward and reverse shock. The observation of the 
dips in different energy bins allows therefore not only the identification of the reverse 
shock but also to trace the propagation of both shocks. Since the neutrino conversion 
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probabilities are energy dependent during the passage of the shocks through the H- 
resonance, neutrino oscillations can be detected even if the energy spectra of different 
neutrino flavors have the same shape but different luminosities. 

We begin in Sec. 2 with an explanation of the reverse shock formation and a 
discussion of our numerical shock-propagation examples. In Sec. 3 we use a schematic 
model of the density profile to derive analytically the generic features imprinted on the 
observable neutrino signal. In Sec. 4 we use a concrete numerical example of shock and 
reverse shock propagation to calculate their typical signatures in a large Ue detector. In 
Sec. 5 we discuss and summarize our findings. 

2. Our numerical SN model 

Relying on the viability of the neutrino-heating mechanism we have performed 
simulations of supernova explosions that followed the evolution in spherical symmetry 
until more than 25 seconds after shock formation. In order to include the effects 
of convection we have also done runs in two dimensions for about one second of 
post bounce evolution. Because of the persistent problems of complete models in 
obtaining successful explosions (see, e.g., Ref. [B]), we replaced the high-density interior 
of the contracting nascent neutron star by a time-dependent inner boundary where the 
neutrino luminosities and spectra were imposed such that the neutrino energy transfer 
to the shock was sufficiently strong for shock revival (for more details, see Refs. [TK|ll6j). 

The density profiles used for the evolving stellar background in the present work (cf. 
Fig. were taken from simulations of a 15 Mq progenitor Although the employed 
neutrino parameters are in the ballpark of results from elaborate neutrino transport 
calculations in supernovae, detailed information about neutrino fluxes and spectra that 
are consistent with explosion models and with the cooling history of the newly formed 
neutron star are currently not available. For this reason we consider the density profiles 
as exemplary for realistic supernova conditions and combine them with a schematic 
description of the neutrino emission that is supposed to contain generic features of the 
expected neutrino signal. 

The space-time evolution of the explosion of this 15 Mq star in a spherically 
symmetric simulation (i.e., neglecting the effects of hydrodynamic instabilities) can 
be seen in Fig. |21 The shock expands from the point of its formation (not shown in 
this figure that displays only times later than 70 ms after core bounce) to a radius 
of about 100 km where it converts into a standing accretion shock. Its subsequent 
slow outward motion reflects the quasi-stationary adjustment of its position in response 
to the decreasing rate of mass infall from the outer stellar layers and to the energy 
deposition by neutrinos in the "gain region" behind the shock. At a radius of about 
500 km the shock reaches a temporary maximum and starts retreating again because of 
decaying neutrino luminosities and a corresponding drop of the neutrino-heating rate. 
The subsequent sharp increase of the neutrino luminosities and of the energy deposition 
behind the shock is sufficient to trigger a successful explosion "in a second attempt." 
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Figure 2. Early evolution (for times between 0.07 and 20 seconds after shock 
formation) of the supernova explosion (with an energy of about 1.5 x 10^^ erg) of 
a 15 Mq progenitor according to a spherically symmetric simulation in which the 
contracting neutron star was replaced by an inner boundary that acts as a neutrino 
source sufficiently intense to cause a neutrino-driven explosion. The black lines 
mark the space-time trajectories of selected mass shells, the blue lines represent 
corresponding contours of constant density, and the thick red line indicates the position 
of the supernova shock. The reverse shock and the contact discontinuity can be 
recognized from kinks and a compression of the isodensity curves in the region behind 
the forward shock. 

This behavior is typical of one-dimensional simulations with conditions close to the 
threshold for shock revival. 

On its way out the shock reverses the infall of the swept-up matter. Continuous 
neutrino energy transfer starts an outward acceleration of heated material in the gain 
layer around the neutron star. At the interface between this "neutrino-driven wind" to 
the outer ejecta a "contact discontinuity" is formed. It is visible as a strip of compressed 
density contours somewhat behind the outgoing supernova shock in Fig. El Even farther 
behind the forward shock, the neutrino-driven wind, whose velocity increases rapidly 
with distance from the neutron star, collides with more slowly moving material and is 
decelerated again. The strongly negative velocity gradient at this location steepens into 
a reverse shock when the wind velocity begins to exceed the local sound speed 'ISj. In 
the model of Fig. Elthis happens at t > 800 ms post bounce at a radius r > 1000 km. 

The forward and reverse shocks are sharp discontinuities where density, pressure (as 
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Figure 3. Radial velocity (upper half panels) and logarithm of density (lower 
half panels) color coded for four different two-dimensional simulations at one second 
after supernova shock formation. The explosion energies of these models are about 
0.3 X 10^1 erg, 0.4 x 10^^ erg, 0.9 x 10^^ erg, and 1.2 x 10^^ erg (from top left to bottom 
right), respectively. The postshock structure is much more complex than in spherical 
symmetry, and depends strongly on the direction. While the contact discontinuity is 
distorted by hydrodynamic instabilities and marks the boundary between low-density 
neutrino-heated bubbles and matter that has been accreted and accelerated by the 
supernova shock, the reverse shock is a clear signature of the onset of a supersonically 
expanding neutrino-driven wind which runs into the earlier, slower supernova ejecta. 
But also multiple (reverse) shocks can occur due to collisions between downflows and 
neutrino-heated, expanding plasma. 



all other state variables) and velocity of the stellar plasma change on the microscopic 
(sub-millimeter) scale of the ion mean free path, on which the dissipation of kinetic 
to thermal energy in the shock occurs. In contrast, the "contact discontinuity" is 
characterized by a density jump but constant velocity and pressure. It emerges in 
the matter in the transition region between shock-accelerated and neutrino-heated 
supernova ejecta. Its width is therefore a result of the explosion dynamics and found 
to be 200-300 km in our one-dimensional model. However, because of our exponentially 
coarsening grid we are able to resolve this outward moving feature only for 2-3 seconds 
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and are unable to make predictions of its exact structure at later times. A schematic 
representation of the described structure is given in Fig. |3J 

The multi-dimensional situation is much more complex (cf. Fig. 12)). The forward 
shock can now be highly deformed, and the reverse shock, which still separates the 
supersonically expanding neutrino-driven wind from subsonically moving outer ejecta, 
shows large angular variation. This reflects the anisotropy of the fluid flow between 
the shock and the neutron star, which is characterized by narrow downflows of cool 
gas and rising bubbles of low-density, neutrino-heated matter. These convective mass 
motions create a highly distorted structure in which multiple shocks can occur due to 
collisions of downflows with expanding, hot matter. Contact discontinuities follow the 
interface between the rising hot bubbles and shock accelerated plasma that has lower 
entropy and higher density. Their shape is affected by Rayleigh- Taylor and Kelvin- 
Helmholtz instabilities. Since the explosion is highly asymmetric, the exact structure 
varies strongly with the direction. Although the generic features of the one- dimensional 
situation are retained, the details of the multi-dimensional structure evolve chaotically 
from random initial perturbations and cannot be uniquely determined by models. 

3. Effect of the shock waves on neutrino propagation 

In present and planned water Cherenkov and scintillation detectors the main neutrino 
detection channel is the inverse beta decay reaction, UeP — >■ ne~^, that allows also for the 
reconstruction of the neutrino energies. Therefore we consider only the z/g spectrum in 
the following. However an analogous analysis can be easily performed in the neutrino 
channel for a detector able to measure the u^. spectrum, using e.g. liquid argon. 

The antineutrino spectra arriving at the Earth are determined by the primary 
antineutrino spectra as well as the neutrino mixing scenario, 

F,^{E) = p{E)Fl{E) + [I - p{E)]Fl{E) , (1) 

where p{E) is the survival probability of a Ve with energy E after propagation through 
the SN mantle:}:, the superscript zero denotes the primary fluxes, and stands for 
either or v^. In general, all quantities of Eq. are not only energy but also time 
dependent. In particular the survival probability p{E) becomes time dependent during 
the passage of the shocks through the H-resonance. Before we discuss the effect of 
the shock propagation on the neutrino spectra, we first recall the case of three-flavor 
neutrino oscillations in a static density profile. In this case, p{E) is not only constant in 
time, but also independent of energy if sin^-i^is < 10~^ or sin^^^is > 10~^. Since about 
60% of the observed neutrinos will be detected in the first two seconds after bounce, 
i.e. before the shock reaches the H-resonance of low-energy neutrinos, it is for many 
purposes sufficient to perform the simpler analysis of neutrino propagation through the 

I We neglect generally Earth matter effects. They slightly increase the chances to detect the SN shock 
propagation, cf. figure 8. 
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sin^ i^is p 
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A 


Normal 


> 10-' 


COS^ 


B 


Inverted 


> 10-3 sin^ 





C 


Any 


< 10-5 sin^-i^Q 


COS^ {}q 



Table 1. Survival probabilities for neutrinos, p, and antineutrinos, p, in various 
mixing scenarios for a static density profile. 

static envelope of the progenitor. For us, the static case discussed in the next subsection 
serves as a reference to measure the modulations induced by the propagating shock. 

3.1. Static case before the arrival of the shock wave at the H-resonance 

The survival probability is determined by the flavor conversions that take place in 
the resonance layers, where pres ~ mj^Am'^ cos2'd/{2\/2GFYeE). Here Am^ and 
are the relevant mass difference and mixing angle of the neutrinos, m-N is the 
nucleon mass, Gp the Fermi constant and the electron fraction. In contrast to the 
solar case, SN neutrinos must pass through two resonance layers: the H-resonance 
layer at pn ~ 10^ g/cm^ corresponds to Am^^^^, whereas the L-resonance layer at 
Pl ~ 10 g/cm^ corresponds to Am^. This hierarchy of the resonance densities, along 
with their relatively small widths, allows the transitions in the two resonance layers to 
be considered independently [Tllj . 

The neutrino survival probabilities can be characterized by the degree of 
adiabaticity of the resonances traversed, which are directly connected to the neutrino 
mixing scheme. In particular, whereas the L-resonance is always adiabatic and appears 
only in the neutrino channel, the adiabaticity of the H-resonance depends on the value 
of '^13, and the resonance shows up in the neutrino or antineutrino channel for a normal 
or inverted mass hierarchy respectively. Table ^ shows the survival probabilities for 
neutrinos, p, and antineutrinos, p, in various mixing scenarios. For intermediate values 
of T^ia, i.e. 10~5 < sin^t^is < 10"^^ the survival probabilities depend even without shock 
on energy as well as the details of the density profile of the SN. 

For large values of {^is and the static density profile of the progenitor, the H- 
resonance is adiabatic. In the case of normal mass hierarchy, scenario A, the resonance 
takes place in the neutrino channel, and antineutrinos are not affected. For an inverted 
mass hierarchy, scenario B, the resonance occurs in the antineutrino channel, so that 
practically all the primary z/g are converted to z/3 and arrive at the Earth as u^- This 
corresponds to an almost complete interchange of Ue and Ux spectra. In scenario C, the 
resonance is strongly non-adiabatic, and hence ineffective. In both the scenarios A and 
C, the primary leave the star as ui, which implies a partial mixing between i>e and 
spectra. 
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Figure 4. Schematic density profile in the presence of a forward and reverse shock 
wave. Between forward and reverse shock a contact discontinuity can appear. 



3.2. Shock waves passing through the H-resonance 

In our numerical model, the outgoing shock wave reaches the H-resonance layer around 
two seconds after bounce and the L-resonance layer around ten seconds after bounce. 
Since for t > 11 s the H-resonance is adiabatic, neutrinos created as arrive in case 
B as P3 at the L-resonance. Thus, the L-resonance that mixes 12-states does not affect 
the survival probability p, independent of the passage of shock waves. In contrast, 
it depends on the neutrino mixing scenario, whether the shock waves modulate the 
neutrino spectra during their passage through the H-resonance: The cases A and C will 
not show any evidence of shock wave propagation in the observed Ue spectrum, either 
because there is no resonance in the antineutrino channel as in scenario A, or because 
the resonance is always strongly non-adiabatic as in scenario C. However, in scenario B, 
the sudden change in the density breaks the adiabaticity of the resonance, leading to 
observable consequences in the Ue spectrum. Let us discuss this scenario henceforth. 

In order to analytically understand the signatures of the shock wave in the observed 
Ue spectrum, let us first assume that the shock wave causes the adiabatic H-resonance 
to become completely non-adiabatic at the forward shock front and the reverse shock 
front, while at all other places the resonance remains adiabatic. In the snapshot of the 
shock wave shown in Fig. I5 non-adiabatic transitions take place only at the two radii 
ri and r2, where the forward and the reverse shock waves are present at that time. 

Neutrinos with energy E undergo H-resonance in the region of density 

The maximum and minimum densities at the two shock fronts, as illustrated in Fig. 
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Figure 5. Survival probability p as a function of energy E assuming \da — i^tl 
(dotted line) or using Eq. Q (solid line) with cos^ = 0.7. 



n/2 



then directly give the range of energies that undergo a non-adiabatic resonance at these 
shock fronts. We introduce now the two approximations of our toy-model used in this 
subsection: First, we assume that all resonances are either completely adiabatic or 
completely non-adiabatic. Second, we require that the density jump at the resonances 
is large, i.e. that the change of the medium mixing angle is close to 7r/2. Then the 
survival probability p{E) is a rectangular wave as a function of energy, as shown in 
Fig. El with a dotted line. Clearly, as the shock wave propagates to lower densities, 
the "turning points" of the square wave shift upwards in energy. The amplitude of the 
square wave, cos^-;?©, is fixed by the value of the solar mixing angle. We can relax the 
second approximation and use the result that, if the density changes sharply from to 
Pa, the survival probability for z/g is 



where and -^b are the values of -^i^ in matter at the densities pa and pb, respectively. 
Using this, the rectangular wave for p gets modified as shown in Fig. El with a solid line. 
However, in order to get a simple understanding of the new features of the conversion 
probability we will use the rectangular wave for p in the rest of this section. 
The primary spectra of neutrinos can be parameterized by PHI 120] 



where (Ei) denotes their average energy, (3i is a dimensionless parameter that relates to 
the width of the spectrum and typically takes on values 3.5-6, depending on the flavor 
and the phase of neutrino emission. Let us approximate the energy dependence of the 
cross section of the inverse beta decay reaction by E"^, and assume an ideal detector. 
Using Eq. (pQ) , Table Q and the survival probability as shown in Fig. the number of 
events observed in a certain time interval is 



p = sin^('(9a — i!}b) cos^ '^q , 



(3) 





A^obs 



Af $ 



PI r(/5pj 



+ cosH^{^,^gl^-^,^gl-) 



(5) 



Signatures of supernova shock propagation 



11 



where A/" is a normalization constant, T{z) is the Euler gamma function, and we have 
introduced the functions 



9k 



Here, T{a,b,c) 



{E.J' 

Ibdt r-i 



El 



(6) 



is a generahzed incomplete Gamma function, and 



a = f3 + k. The values Ei, E2, E^, E4 are the same as Eia, E^, E2a, E2b in ascending 
order. The mth moment of the total energy E^^ deposited in the detector during a 
certain time interval is similarly given by 



ZT"' 



obs 



AT 



{E,J^+^T{(3,^+2 + m) 



2+m 



+ COs2 790($p^(72';m-^P.^?2;. 



(7) 



Since the values of Eia, Eib, E2a, E2b change as a function of time, so do the values 
of Nohs,Eobs- As time goes on, the shock wave reaches layers with smaller densities 
and therefore, according to Eq. higher-energetic neutrinos are affected by the wave 
propagation. The time dependence of the numerical profiles can be imitated starting 
from t = Is by Ei{t) = (t/s)"' \n{t/s) MeV, and where = {1.3,1.5,1.7,2.0}. As a 
consequence, the double square structure in the survival probability shown in Fig. 
keeps on shifting towards higher energies with time. 

In order to study the model dependence of our results, we consider three models 
that give very different predictions for the neutrino spectra. Two of them, Gl and 
G2, are motivated by the recent Garching calculation |2I] that includes all relevant 
neutrino interaction effects like nucleon bremsstrahlung, neutrino pair processes, weak 
magnetism, nucleon recoils and nuclear correlation effects. The third one is from 
the Livermore simulation [22] that represents more traditional predictions for flavor- 
dependent SN neutrino spectra used in many previous analyses. The parameters of 
these models are shown in Table |21 



Model 






{Eo{>^x)) 






L 


12 


15 


24 


2.0 


1.6 


Gl 


12 


15 


18 


0.8 


0.8 


G2 


12 


15 


15 


0.5 


0.5 



Table 2. The parameters of the used primary neutrino spectra models motivated from 
SN simulations of the Garching (Gl, G2) and the Livermore (L) group. We assume 
— 4 and /3p^ = 5. 



In Fig.ini we show the values of the average event energy (Ee) = -Eobs/A'obs observed 
(left) as well as ^ = {E"^) / {E^)"^ (right) as functions of time for the three models of Tab|^ 
The parameter ^ measures the pinching of the observed spectrum and indicates therefore 
if it is a superposition of primary spectra with different shapes. We compare the case 
where both forward and reverse shock waves are present with the forward shock only 
case. We approximate the latter by setting E2a = -£'26- 

In the absence of shock waves, the adiabaticity of the propagation in case B ensures 
that the z/g flux arriving at the Earth, Fp^{E), basically coincides with Fp^{E). Thus, 



Signatures of supernova shock propagation 



12 



breaking that adiabaticity will imply a partial presence of F^^{E) in Fp^{E). Therefore 
any departure of the observables from their values in the no-shock case induced by the 
shock depends strongly on which part of the energy spectra is affected by the shock 
wave, and how different the original spectra are at those energies. 

At early times, t < 4 s, the shock affects only the low-energy part of the spectra, 
E <20 MeV. At these energies, the ratio [E) / Fp^ (E) varies for different SN models: 
it is smaller than one in model L and larger than one in models Gl and G2. Therefore, 
any modulation in the observables at these early times will strongly depend on the 
model considered. At later times, t > 4 s, neutrinos with higher energies, E >20 MeV, 
feel the breakdown of the adiabaticity. At these energies, the Ux flux is in all models 
larger than the z/g flux, F^^{E) > F^^{E). Therefore our prediction for the observable 
modulation of SN neutrino spectra at later times, t > 4 s, will be model independent 
and thus more robust. For this reason, we will basically concentrate the discussion on 
the modulation of the observables arising at these times. 

The main effect visible in the time evolution of {Ee) shown in the left panel of 
Fig- El is a decrease of {E^.) when the shock passes the resonance region and reduces the 
number of events at high energies. However, the details will depend on the features of 
the shock propagation. If only the forward shock wave is present, a deep and wide dip 
is expected. If moreover a reverse shock is formed, a double dip will be observed: the 
positions of the two dips correspond to the times when one of the peaks in the survival 
probability in Fig. El coincides with the peak energy of the antineutrino spectrum. In the 
case of ^ = {El) / {Ef.)"^ the double dip structure becomes a double bump. In contrast 
to an energy-independent survival probability, ^ is not always increased by mixing but 
can be also reduced. 

It is remarkable that even in the model 02, where the original mean energies {Ei) 
are exactly the same, the energy dependence of the conversion probability results in an 
observable effect of the shock wave propagation, because of the flavor-dependent fluxes. 

Note that in this section we have assumed that the level crossings are completely 
adiabatic except at the forward and reverse shocks. For sin^ t^is close to its current upper 
limit, this is very likely to be the case. In the next section, we numerically calculate 
the survival probabihty without the foregoing assumption, and find that the qualitative 
features obtained here do indeed survive in the complete numerical calculation. These 
features are thus robust, and though illustrated here for the case of a special shock 
profile, should be valid for a general case. 

4. Detection of the shock and the reverse shock 

In this section, we examine in more detail the signatures of the forward and reverse 
shock propagation imprinted on the neutrino signal of a future galactic SN. The different 
ingredients used for the simulation are 

• Initial fluxes: we assume for the time dependence of the total luminosity the one 
given in Ref. [221, shown in the left panel of Fig. [3 We use the parameterization 
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Figure 6. Time dependences of tlie observables (Ee) (left) and {E^)/{Ee)^ (right) for 
different primary spectra, (L, Gl, G2), and assuming both forward and reverse shock 
(blue solid line) and only forward shock wave (red dashed line). 




Figure 7. Left: Time dependence of the neutrino luminosities L used, from Ref. |22| . 
Right: Density profiles at different times obtained from the numerical simulations of 
section 2. The H-resonance density, pn, for different neutrino energies, 5, 30 and 60 
MeV, is also shown. 



of Eq. for the spectral shape, with the values given in Tab. El 

• The density profile of the SN envelope at different times as discussed in Sec. EI 
In the right panel of Fig. [71 we show the effect of the forward and the reverse- 
shock wave on the density profiles at several instances. The basic structures of the 
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numerical profiles coincide well with the schematic profiles from Fig. 0] used in our 
analytical discussion, except at early times, where p2a ^ Pia- 

• The numerically calculated survival probability p{E,t). 

• Detector: we simulate the neutrino signal at a megaton water Cherenkov detector 
assuming a SN distance of 10 kpc. The detector response is taken care of in the 
manner described in Ref. j2S]- The neutrino signal is dominated by the inverse 
beta reaction u^P ne~^, while all other reactions have a negligible influence on 
the analysis below. In future, the addition of Gd may allow for the efficient tagging 
of this reaction [21] . Therefore, we take into account only the inverse beta reactions 
in the following analysis. 

• Neutrino mixing parameters: we concentrate on case B, i.e. an inverted mass 
hierarchy and large t^is. In particular, we use the following numerical values 
for the relevant parameters, Am^ = 6.9 x 10~^ eV^, tan^t?© = 0.42, Aml^^ = 
2.6 X 10~^ eV^, and tan^-i^is = 10~^. However, we have checked that the signatures 
of the shock propagation like the "double-dip" feature are detectable even for 
tan2^9i3 = 5 X 10"^ 

The key ingredient to observe signatures of the shock wave propagation is the time 
and energy dependence of the neutrino survival probability. In Fig. |H1 we show the 
energy dependence of p{E, t) averaged with the energy resolution function of Super- 
Kamiokande, for the case with (right panel) and without (left panel) a reverse shock. 
To simulate the latter case, we removed by hand the reverse shock from the numerical 
density profiles. 

Qualitatively, we observe the same features as using the analytical approximation 
of Sec. El While the survival probability has a single peak if there is only a forward 
shock propagating, the reverse shock superimposes a dip at those energies for which 
the resonance region is passed by both shock waves. All these structures move in time 
towards higher energies, as the shock waves reach regions with lower density. The 
width and location of the peaks depend on the particular values of Eia, Eib, E2a, E2b at 
each instant. We can observe, for example, that at t = 2 s only low energy neutrinos, 
~ 10 MeV, are affected by the reverse shock. A little bit later, at t = 5 s, neutrinos 
with energies E ~ 10-15 MeV are influenced by the forward shock wave (the first peak), 
while those with E ~ 20-60 MeV feel the effect of both the reverse and the forward 
shocks. The effect of the two strongly non-adiabatic transitions cancel partially and 
p{E, t) has a valley around ~ 30 MeV. For > 60 MeV, neutrinos cross only one shock 
wave, thus a second peak appears in p{E,t). At later times, this pattern is basically 
repeated, only shifted to higher and higher energies. 

The first observable we consider is the average of the measured positron energies, 
{Ee), with a time binning of 0.5 s. In Fig. IHl we show (E^) for various cases together 
with the one sigma errors: Each of the six panels contains the case that no shock wave 
influences the neutrino propagation, the case of only a forward shock wave and of both 
forward and reverse shock wave. In the three left panels the (E^) are constant in time, 
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Figure 8. Survival probability t) as function of energy at different times averaging 
in energies with the energy resolution of Super-Kamiokande; for a profile with only a 
forward shock (left) and a profile with forward and reverse shock (right). At i = 5 s, we 
show p{E, t) including Earth matter effects for a zenith angle of the SN of 62° (black 
line). 



while in the three right panels we assume that the (Ei) decrease linearly after 5 s. 
Finally, the upper, middle and lower panels show the three models L, Gl and G2 for 
the primary neutrino spectra given in Table El 

As our most important result we find that the effects of the shock wave propagation 
are clearly visible in the neutrino signal of a megaton detector, independent of the 
assumptions about the initial neutrino spectra. Moreover, it is not only possible to 
detect the shock wave propagation in general, but also to identify the specific imprints 
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Figure 9. The average energy of Dp — *■ ne^ events binned in time for a static density 
profile (magenta), a profile with only a forward shock (red) and with forward and 
reverse shock (blue). In the left panels we use (Ei) = const., whereas in the right 
panels (Ei) decreases for t > 5 s. The error bars represent 1 a errors in any bin. 
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Figure 10. Left: Time dependence of (Ee) for a profile with a forward and reverse 
shock for several values of tan^ as indicated; for model L. Right: Time dependences 
of the observable ^ = {El) / {Ef.)'^ for a static density profile (magenta), a profile with 
only a forward shock (red) and with forward and reverse shock (blue) . 



of the forward and reverse shock versus the forward shock only case. The signature of 
the reverse shock is its double-dip structure compared to the one-dip of a forward shock 
only, as predicted by our analytical model in Fig. IHl While the presence or absence of 
structures at early times does depend on the initial fluxes and the details of the density 
profiles, the two-dip pattern is a robust signature of the presence of a reverse shock. 
To study the dependence of the double-dip structure on the value of t^is, we show {E^.) 
as function of time for different 13-mixing angles in the left panel of Fig. Even for 
as small values as tan^ ^9i3 = 5 x 10~^ the double-dip is still clearly visible, while for 
tan^ di^ = 1 X 10~^ only a bump modulates the neutrino signal. 

Next we discuss the possibility to detect the imprint of the shock waves in other 
observables. To shorten the exposition, we consider only one model (L) together with 
a linear decrease of {Ei) after t = 5 s. In Fig. EH "we show the time dependence of the 
observable ^ = {ED/^Ee)"^. If the shock does not influence the neutrino spectra, i.e. 
p ~ 0, then ^ = ,^0 = + 3)/(/3p^ + 2). For a mixed spectrum, with p{E) independent 
of energy, we have ^ > ^0 iii general. However, for a strongly energy dependent p{E) as 
in the case of a shock, can increase or decrease depending on the details of p{E) and 
the primary spectra. Such modulations in the value of ^ that coincide in time with the 
modulations in {E^,) act as a confirming evidence for the passage of the shock wave. 

Another observable that can be exploited is the total number of detected events. In 
the left panel of Fig. ^2 we show the number of events binned in time. This observable 
is particularly interesting for detectors like IceCube with poor or no energy resolution at 
all for SN neutrinos. If the luminosities are fast decreasing with time, the modulations 
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Figure 11. Left: Total number of events detected as function of time. Riglit: Number 
of events binned per energy decade as function of time for forward and reverse shock. 



introduced by the time-behavior of p{E, t) may be difficult to disentangle from the 
overall time-behavior of the luminosities. Nevertheless, even a detector without energy 
resolution like IceCube has a potential to discover the imprint of propagating shock 
waves on the neutrino signal, if the luminosities are not decreasing as fast as found in 
Ref. 

Since the information about the shock fronts is encoded in the energy and time 
dependence of the survival probability p{E, t), one can expect that the number of events 
in a fixed energy range is a more useful observable than the total number of events. In 
the right panel of Fig. ^2 we show therefore the number of events binned in energy 
intervals of 10 MeV as function of time for the case of a reverse shock. We can observe 
clearly how the positions of the two dips change in each energy bin. It is remarkable 
that the double-dip feature is not only stronger after binning, but also allows one to 
trace the shock propagation: Given the neutrino mixing scheme, the neutrino energy 
fixes the resonance density. Therefore, the progress of the shock fronts can be read off 
from the position of the double-dip in the neutrino spectra of different energy. 

We illustrate this by examining the time evolution of the number of events in the 
bin 40 ± 5 MeV in detail. From Fig. El right panel, it can be read off that between 
4.5-7.5 s, the neutrinos with this energy pass through two nonadiabatic resonances, 
i.e. both the forward and the reverse shock are present in the regions with density 
P40 ~ 850 g/cm^. Between 7.5 s and 9.5 s, only one of the shock fronts is present at this 
density. This inferred behavior of the shock wave can be seen to correspond with the 
shock proffie used, see the left panel of Fig. El For times before 4.5 s and after 9.5 s, 
the data is not able to say anything concrete about the shock wave propagation. 
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Figure 12. Left: Density profiles at different times witli tlie resonance layer of 
neutrino with E = AO MeV superimposed. Right: Number of events expected in the 
energy range = 40 ± 5 MeV as function of time. 



If only the forward shock is present, then Fig. still allows us to infer that one 
shock wave was present in the region p^Q between 5-9.5 s. No concrete conclusions can 
be drawn about the behavior of the shock wave beyond this time interval. 

The number of shock waves present in a region with any given density ~ 10^ g/cm^ 
can similarly be extracted from the data by considering the time evolution of the number 
of events in the energy bin corresponding to that density. An extrapolation would allow 
one to trace the positions of the forward and the reverse shock waves for times between 
1-10 s. 

Until now we have neglected the possibility that other discontinuities in addition to 
the forward and reverse shock can influence the neutrino propagation. However, because 
of the decreasing resolution in our simulation at large radii as described in Sec. 2, we are 
only able to resolve the outward moving contact discontinuity for 2-3 seconds. Its width 
at that time is found to be 200-300 km in our one-dimensional model. Hydrodynamic 
instabilities are expected to widen the contact discontinuity, and we expect therefore 
that the neutrino evolution is adiabatic also at later times at this discontinuity. But 
to get an idea of the possible influence of this discontinuity on observables, we have 
steepened it by hand for t > 3 s, so that the neutrino evolution becomes strongly non- 
adiabatic at this point. In Fig. [T31 we show the detected {E^) in this case. Since now 
there are three shocks (and six associated energies), the structures in {Eg) become more 
complicated but the imprint of the shock wave is still clearly visible. In particular, 
there are effects of the shock wave propagation at later times than in the previous cases, 
because for the density profiles considered in our study the largest density affected by 
this new discontinuity is larger than the maximum of pia and pn,. 
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Figure 13. Left: Observed (Ee) for a static density profile (magenta), a profile with 
only a forward shock (red), with forward and reverse shock (blue), and in the case 
that the contact discontinuity is steepened by hand (green). Right: Observed (Ee) in 
Super-Kamiokande, color code as in the left panel. 



Finally, we want to discuss if the largest currently operating detector, Super- 
Kamiokande, can observe the signatures of shock wave propagation. Compared to a 
megaton detector, the event rate in Super-Kamiokande is reduced by a factor of 30. 
Therefore, we plot in the right panel of Fig. the mean event energy {E^) with one 
second time bins in order to decrease the statistical error. The features of the different 
shock propagation scenarios remain visible, but are now more difficult to disentangle 
because of the larger errors. However, one should keep in mind that the number of 
events is not only related to the volume of the detector but also to the distance to the 
SN. Thus, the potential of Super-Kamiokande to yield important clues about the shock 
propagation cannot be dismissed. 

5. Summary 

We have performed simulations of supernova explosions adjusting the neutrino energy 
transfer to the shock so that a successful explosion resulted. Following the time evolution 
in spherical symmetry until more than 25 seconds after shock formation, we have found 
that in addition to the forward shock wave a reverse shock forms when the supersonically 
expanding neutrino-driven wind collides with the slower earlier SN ejecta. Both the 
forward and reverse shock are sharp discontinuities where the density changes on a 
sub- millimeter scale. 

The sudden density jump at the two shock fronts results in a strongly non-adiabatic 
evolution of neutrino flavor oscillations, when the shocks cross the H-resonance layer. For 
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a "large" 13-mixing angle, sin^'!?i3 ^ 10~^, the MSW enhancement of flavor conversion 
that is otherwise working is interrupted during the shock passage. This break-down 
of adiabaticity results in a reduction of the average energy and number of detected Ue 
events for an inverse neutrino mass hierarchy (case B) or of events for a normal 
hierarchy (case A). The reduced event rate during the shock passage through the H- 
resonance allows even experiments with missing energy resolution for SN neutrinos, as 
e.g. IceCube, to observe effects of the shock propagation. 

The characteristic signature for the presence of two shocks is the "double-dip" 
feature in the time-binned energy spectrum of observed electron-neutrino events. We 
have found that already Super-Kamiokande can potentially distinguish between different 
shock propagation scenarios, while a megaton detector may be able to trace the time- 
evolution of forward and reverse shocks. 

The modulation of SN neutrino spectra by propagating shocks is only possible 
if t?i3 is "large," sin^i?i3 3> 10~^. Observing these features in the neutrino signal of 
a future galactic SN would therefore indicate that the measurement of ■J?i3 and the 
detection of leptonic CP violation may be within reach. In addition, the detection 
of any modulation in the Pe or the Ue spectrum by shock effects would identify the 
neutrino mass hierarchy. Remarkably, observing features of SN shock propagation is 
complementary to the detection of Earth matter effects: while observing modulations 
by SN shocks in the Pe spectrum identifies case B, modulations by Earth matter effects 
exclude this case. 

We want to close with a speculative remark about the time-structure of the neutrino 
events from SN 1987A. The events observed by the Kamioka experiment can be grouped 
into two time bins: nine events during the first two seconds, then three more events 
after a time gap of six seconds. Intriguingly, such a time structure is compatible with 
the modulation of the neutrino signal by SN shocks, if the neutrino luminosities are 
decreasing slowly with time. Thus, one might speculate that this gap is connected with 
the passage of SN shocks through the H-resonance. Then, the SN 1987A signal would 
be a hint that case B is the true neutrino mixing scenario, i.e. that the neutrino mass 
hierarchy is inverted. 
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